Ultrasonic attenuation in d-wave superconductors 
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The longitudinal ultrasonic attenuation at low temperatures for a ei-wave superconducting order pa- 
rameter is examined. The ratio of the superconducting to normal state attenuation, a 3 /a n , depends 
strongly on the direction of propagation of ultrasound with maxima corresponding to nodes of the 
order parameter on the Fermi surface. We propose that measurements of ultrasonic attenuation in a 
single crystal of the high-T c superconductors can give additional information regarding the pairing 
symmetry. 

PACS numbers: 74.20.-z, 74.40.+k, 74.50. +r, 74.70.Mq 

An important aspect of the high-T c superconductors 
that is currently receiving a great deal of attention is 
the angular dependence of the energy gap. Some ex- 
perimental evidence [[l]j|] has been brought forward for 
d-wave pairing instead of the more conventional BCS s- 
wave pairing. The purpose of this note is to point out 
that ultrasonic attenuation may offer a very effective and 
useful tool for the experimental determination of the vari- 
ation of the energy gap around the Fermi surface. We 
demonstrate that, due to the layered structure of the 
high-T c cuprates, the anisotropy of the d-wave order pa- 
rameter should result in a very pronounced dependence 
of the ultrasonic attenuation on the direction of propaga- 
tion (parallel to the C11O2 layers) at low temperatures. 

In order to concentrate on our main point, we as- 
sume that the scale of the quantum of sound energy, 
Tiuj qi is negligibly small compared to the energy gap. 
Typical experimental values of uj q are of the order of 
10 9 — 10 10 sec" 1 w 0.01 — 0.1 meV, equivalent to a tem- 
perature of 1 K or less. (For simplicity, we ignore the im- 
mediate vicinity of a node.) This assumption leads to the 
picture of the ultrasonic attenuation as a kind of quasi- 
elastic scattering of the quasiparticles of energy E by the 
incident sound quantum. The initial states are counted 
in terms of increments of the normal state energy, de = 
(de/dE)dE. The integration for the total transition rate 
in the Fermi "golden rule" formula needs to include an 
additional factor de/dE for the final state density. But 
these two factors of de/dE are cancelled by the squared 
strength of the coupling of the sonic field with the quasi- 
particle energy j|, which is proportional to dE/de. This 
follows from the fact that it is to e, and only indirectly to 
E, that the lattice deformation is coupled. Consequently, 
the ratio of the longitudinal ultrasonic attenuation in the 
superconducting and normal state is given by the formula 
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X = v~'cos 26 
FIG. 1. The scaling function for the angular dependence of 
the d-wave order parameter. This function characterizes the 
universal behavior which connects small doping and the angle 
dependence of the order parameter. 



sound quantum is not negligibly small, as can be seen in 
the standard treatment [0. When A differs, as presently 
under consideration, for different quasi-particles at the 
Fermi surface we take the average of the right hand side 
of Eq. (1), weighted by the contribution in the normal 
state. It is possible to obtain the same result using 
Green's functions to calculate the density-density correla- 
tion function, the imaginary part of which is proportional 
to ultrasonic attenuation. The result for longitudinal ul- 
trasonic attenuation in anisotropical superconductors is 
given by the following expression Q : 



a s (0) < <S(v k ■ q) |.g(k,k)| 2 2/(A k ) > FS 



a n (9) 



<<S(v k -q) |5(k,k)| 2 > FS 



(2) 



df(E) 
dE 



dE = 2/(A) 



(1) 



exp (3 A + 1 

where f(E) is the Fermi function and A is the supercon- 
ducting energy gap. This derivation || of Eq. (1) does 
not involve the BCS coherence factors. They do, how- 
ever, enter for frequencies so high that the energy of the 



where g(k, k) is the matrix element of electron- (acoustic) 
phonon interaction. The momentum k and velocity v k 
are on the Fermi surface and < ... >fs is the integration 
over the Fermi surface. The angle 9 describes the direc- 
tion of q, the transfer momentum. In the case of a two- 
dimensional Fermi surface, the integrations in Eq. (2) 
are reduced because of the (5-functions and the matrix 
element g(k, k) in Eq. (2) cancells. For concreteness, we 
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FIG. 2. The dependence of longitudinal ultrasonic atten- 
uation at T/T c — 0.1 and different band-filling factors. 
The dotted line shows the analytic result at n = 0.45 (see 
Eq. (5-6)). 



consider the two-dimensional tight-binding electron spec- 
trum 



£k = — 2t [cos k x + cos k y ] — fi . 



(3) 



where t is the hopping matrix element, [i is the shift 
of chemical potential as determined by the band-filling 
factor n (the Fermi surface is specified by ek = 0). The 
lattice constant is taken to be 1 for simplicity. The d- 
wave order parameter can be written in the form 



Ao(T) 



cos k x — cos k y ] 



(4) 



where Ao(T) obeys the usual BCS temperature depen- 
dence. One can see that, in any given direction, the ultra- 
sonic attenuation falls off exponentially with temperature 
except for directions close to those of the nodes of the d- 
wave order parameter. Exactly at the nodes, a s /a n = 1, 
independent of T. This yields sharp spikes in the angular 
dependence for a s /a n , which become more pronounced 
with decreasing temperature. The analytic result for the 
case where the band-filling factor is close to one-half ( 
n ~ 0.5 - 3ir- 2 v[l + ln(2/^)], with v = |/i|/4f < 1) is 
a kind of universal scaling, so that all such small values 
of v are described by the single scaling function 



F(X) 



X 



(5) 



which is plotted in Fig. |l| versus the scaling variable \ — 
f -1 |cos20|. In this case of the square lattice, 9 is the 
angle between q and the y-axis. Thus, Eq. (2) becomes 



s(0) 



a n (6) 



1 + CXP (M^ W ) 



(0) 



It follows from Eqs. (5) and (6) that, in the low- 
temperature region (T <C T c ), very sharp peaks in 
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FIG. 3. The angle dependence of the ratio of longitudi- 
nal ultrasonic attenuation in the superconducting state to the 
normal state value at n = 0.4 and different temperatures. 
Note the exponential decrease of ultrasonic attenuation with 
temperature far from nodes, as well as the closeness to one of 
the ratio a s /a n near the nodes. 



a a (9)/a n (6) appear at 6 = ±45°, ±135°. The sharp- 
ness of the peaks is the result of the anisotropy of both 
the gap and the Fermi surface. The dotted curve in Fig. 
H illustrates the application of Eqs. (5) and (6) at the 
low temperature T/T c = 0.1 for n = 0.45 (or v = 0.065). 
The solid curve is calculated numerically for the same fill- 
ing without recourse to the small v approximation. The 
closeness of the solid and dotted curves is an indication of 
the accuracy for v = 0.065, of this approximation and of 
the scaling function that follows from it. The numerical 
computations for other fillings are exhibited by the other 
three curves in Fig. ^. These curves clearly demonstrate 
the salient feature mentioned above, namely, that at low 
temperature, and even at fillings significantly far from 
half- filling, a sharp peak can be expected in a s /a n at 
the sonic propagation direction that corresponds to the 
nodes in the energy gap. Figure [| shows the development 
of this peak, for n — 0.4, as the temperature is lowered. 

To conclude, we note that effect of anisotropy of the or- 
der parameter upon ultrasonic attenuation was observed 
in various metals || and in superfluid 3 He [[r). Some re- 
cent progress has been made in ultrasonic measurements 
of single crystals of UPtj, || , where a three-dimensional 
d-wave order parameter is assumed. The earliest ultra- 
sonic measurements on high-T c crystals found a quite 
small velocity change at the superconducting transition 
temperature, which indicates that the influence of super- 
conductivity on the acoustic processes is not yet fully 
understood ||. Nevertheless, the sharp spike- like nature 
of the angular dependence of the longitudinal ultrasonic 
attenuation in the case of the possible d-wave order pa- 
rameter should be observable and should provide a clear 
signature for the latter. 
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